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QQ phase a BCS superfluid of diquark pairs is formed, and the pattern of chiral symmetry breaking 



is very different from that at low density. Analytical results for the spectral density obtained from 
this random matrix theory allow for the extraction of the BCS gap from lattice data. 
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1. Introduction 

Lattice studies of QCD at nonzero quark chemical potential /i are hindered by the infamous 
sign problem, see [1] for a review. Two-color QCD with an even number of pairwise degenerate 
quarks does not have a sign problem and can therefore be simulated on the lattice [2] . It shares 
many qualitative features, such as confinement and chiral symmetry breaking, with three-color 
QCD, but the detailed features of both theories are rather different, such as the pattern of chiral 
symmetry breaking, the particle spectrum, or the phase diagram. Nevertheless, two-color QCD is 
an interesting theory in its own right. It has been studied in great detail at zero and low density, 
see, e.g., [3]. In this contribution, we also address the region of high density in which the pattern 
of chiral symmetry breaking is different from that at low density and in which a BCS superfluid of 
diquark pairs is expected to be formed because there is an attractive channel between quarks near 
the Fermi surface. In earlier work [4] , we have derived the low-energy effective chiral Lagrangian 
for jj. » Asu(2)> identified the corresponding £-regime, and derived Leutwyler-Smilga-type sum 
rules for the eigenvalues of the Dirac operator. This work has been summarized at Lattice 2009 [5]. 

In the lowest order of the £-regime, sometimes also called "microscopic domain", the theory 
becomes zero-dimensional. This zero-dimensional limit of the theory can alternatively be described 
by a random matrix theory (RMT). Many examples of such an exact mapping are known, in par- 
ticular for two- and three-color QCD at zero and low density, see [6, 7] for reviews. Therefore 
the natural question is what random matrix theory describes the microscopic domain of two-color 
QCD at high density. The answer to this question was given in [8] and will be reviewed in Sec. 2. 
The advantage of having a random matrix theory is that it allows us to compute a large number of 
analytical results characterizing the Dirac eigenvalues, see Sec. 3. This task would be much more 
difficult in the effective theory. The analytical results at high density contain the BCS gap A, which 
was computed for asymptotically high density in a weak-coupling approach in [9, 10], as a param- 
eter. Therefore, A can be extracted from lattice data for the Dirac eigenvalues by matching them to 
the analytical results from random matrix theory. Another interesting feature of two-color QCD is 
that it allows us to study the sign problem, either for an odd number Nf of flavors by turning on /i, 
or for even Nf by detuning the quark masses from their degenerate values, see Sec. 4. 

2. Random matrix theory at low and high density 

The pertinent random matrix theory for two-color QCD at low density has been formulated in 
[1 1], with partition function 

Zrmt(m) = / dPdQe-^MP^'^QQ'^Yi^ct ( J'l ^^^^\ , (2.1) 
J fj^ y-P' +nQ' nif j 

where the ntf are the quark masses, P and Q are real matrices of dimension N x [N + v), dP and 
dQ we. Cartesian integration measures, N is assumed to be proportional to the Euclidean space- 
time volume V4, and v can be identified with the topological charge. Note that sometimes other 
conventions for the width of the Gaussian distribution of P and Q are used in the literature. The 
RMT Dirac operator D{}x) is the matrix in (2.1) with the mass set to zero. 
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It was shown in [8] that in the N ^ oo Umit this RMT partition function is identical to the par- 
tition function obtained from the static (or zero-dimensional) effective Lagrangian for two-color 
QCD at low density^ given in [3]. More precisely, the two partition functions have the same de- 
pendence on the quark masses and on the chemical potential. The mapping between dimensionless 
RMT quantities and physical quantities is given by 

^/Nm = mphysGV4 = m^F^Vi , (2.2a) 



1 



-^M =i"phys^'V4, (2.2b) 

where G and F are low-energy constants in the effective Lagrangian in the notation of [3]. 

Both the random matrix theory (2.1) and the con^esponding effective Lagrangian explicitly 
depend on the chemical potential /i. In contrast, the effective Lagrangian at high density derived 
in [4] does not explicitly depend on jj.. It only depends on the quark masses, which appear in the 
combination m^A^V4. A hint as to what the correct random matrix theory at high density should be 
can be obtained by noting that (2.1) is basically symmetric under /i — 1/jLt (except that real and 
imaginary parts are interchanged). Maximum non-Hermiticity, which is expected at high density, 
corresponds to jLt = 1 . We therefore conjecture that at high density the random matrix theory is, 
after a redefinition of the random matrices, given by 



Zrmt = / ftdet (^^/ , (2.3) 



where the dimension of A and B is again N x {N + v)? In the high-density phase, we restrict 
ourselves to an even number of flavors. 

Let us first check that we obtain the correct pattern of chiral symmetry breaking. To this end, 
we rewrite the A/y-flavor determinant resulting from (2.3) in the chiral limit in the form 



det^/fVVdet^^/^f '/Vet^'/ 




(2.4) 



The matrices in the two factors on the RHS of this equation have the form of the chiral orthogonal 
ensemble of random matrix theory. It was shown in [12] that the symmetry breaking pattern in that 
ensemble with Nf/2 flavors is V{Nf ) — )• Sp{Nf). Since we have two such factors, (2.3) with Nf 
flavors has the symmetry breaking pattern \J{Nf) x \J{Nf) — )■ Sp{Nf) x Sp(A'y ). This agrees with 
the symmetry breaking pattern in the effective theory due to the formation of a diquark condensate, 
which is given by SU(A^/)l x SU(A^/)r x U(1)b x U(1)a ^ Sp(A^/)l x Sp(A^/)i? [4]. 

We have also shown [8] that in the N ^ oo limit the RMT partition function (2.3) is identical 
to the partition function of the high-density effective theory in the zero-dimensional limit, i.e., the 
two partition functions have the same mass dependence. The mapping between the dimensionless 
RMT mass and the physical mass is now quite different from (2.2a) and given by 

V3 r- 

m = — mphysAy V4 . (2.5) 

71 



^By low density we here mean the regime of weak non-Hermiticity, see Sec. 3 for the definition of this regime. 
^Only the case v = is physically relevant since topology is strongly suppressed at high density. 
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The arguments presented so far, while giving overwhelming evidence in favor of the equiva- 
lence of the random matrix theory (2.3) and the effective theory at high density, do not constitute 
a full proof. For such a proof one would have to show that all spectral coiTclation functions are 
identical in both theories, which requires studying the partially quenched version of the theory. 
Such a study has not been done yet, but we have no doubt that the outcome would be positive. 

3. Exact results from random matrix theory 

We can now proceed to compute spectral correlation functions from the random matrix theory 
in the N ^ °° limit. At /i = the RMT eigenvalues A are purely imaginary, while at /i 7^ 
they are either purely real, purely imaginary, or come in complex conjugate pairs. We are mainly 
interested in the so-called microscopic spectral density of the small eigenvalues, i.e., we rescale all 
eigenvalues by a quantity 5 that is, up to a numerical prefactor, equal to the mean level spacing 
near zero. This results in complex numbers z = A/5 of order 0(1). To see an effect of the quark 
masses on the small eigenvalues we need to rescale them in the same way, resulting in my = my /5. 

The random matrix theory can be solved in two different regimes: 

• In the regime of weak non-Hermiticity, the combination fi^ = INpL^ = Afi^^^^F'^'V^ is kept 
fixed in the limit N ^ °°. While this regime might appear to be mainly of academic interest 
since /X — )• in the thermodynamic limit, it has an important phenomenological application, 
i.e., the extraction of the low-energy constants G and F from lattice data. In our conventions 
we have 5 = 1 jl^f^ in this regime, i.e., z = 2\/NX and m = 2y/Nm. 

• In the regime of strong non-Hermiticity, pL is kept nonzero in the limit N The analytical 
results in this regime are the /t — 00 limits of the corresponding weak non-Hermiticity results. 
Their functional form is identical for all < /X < 1, and the jLt -dependence only enters through 
a rescaling of the eigenvalues. In our conventions we have 5 = 1 in this regime so that no 
A^-dependent rescaling of the eigenvalues and the masses is necessary. 

In [13] analytical results for the microscopic spectral "density" (which we put in quotation 
marks since this quantity can become negative if there is a sign problem) were obtained in both of 
the above-mentioned regimes in the quenched case, i.e., for Nf = flavors. In the meantime, the 
generalization to the unquenched case has been worked out [14]. Since the analytical results are 
rather cumbersome we will not present them here. Important features of the results are exhibited 
in Figs. 1 through 4. Comments on these features are given in the figure captions. 

4. The sign problem 

As a good measure of the sign problem in two-color QCD we define the quantity 




(4.1) 



see [14] for a more detailed discussion. 
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Figure 1: Microscopic spectral "density" of the purely imaginary (left) and purely real (right) eigenvalues 
in the regime of weak non-Hermiticity foT Nf = I, fL — 3, V — 0, and different values of m. The density of 
the real eigenvalues goes through zero for x = nif. 
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Figure 2: Microscopic spectral "density" of the complex eigenvalues in the regime of weak non-Hermiticity 
for Nf — 1 and v = 0. Left: /i = 1.8 and m ~ 0. The massless quark causes a depletion of the density 
near the origin. Right: fx ~6 and m — 20. For large fi there is an elliptical domain in which the "density" 
oscillates strongly. 
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Figure 3: Microscopic spectral "density" of the purely imaginary (solid) and purely real (dashed) eigenval- 
ues in the regime of strong non-Hermiticity for Nf ~ 2, m\ ~ 2, 1112 = 3, V = (left), and v — 2 (right). The 
density of the real eigenvalues goes through zero for x = mf. 




Figure 4: Microscopic spectral "density" of the complex eigenvalues in the regime of strong non-Hermiticity 
for Nf = 2 and v = 0. Left: mi —1712 — 2. Degenerate masses generate a dip in the spectrum at z = m. Right: 
mi = 2, 1112 = 8. Unequal masses result in a domain of strong oscillations, indicative of the sign problem. 
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Figure 5: Average sign at weak non-Hermiticity for A^/ = 1 as a function of (1. Left: v = is kept fixed and 
m is varied. Right: m = is kept fixed and v is varied. 
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Figure 6: Average sign at strong non-Hermiticity for Nf — 2, thi — I, and v = 0, 1 , 2 as a function of m2- 

We first consider tlie regime of weak non-Hermiticity. We clioose Nf = 1 and turn on fi to 
study its effect on the average sign, see Fig. 5. It is evident from the plots that the sign problem 
(i) increases with increasing jl, (ii) decreases with increasing m, and (iii) decreases with increasing 
V (in agreement with [15]). A quantitative analysis [14] reveals that in the thermodynamic limit 
the average sign makes a first-order transition from 1 to at /t = which in physical units 

corresponds to a critical chemical potential /Iphys = mji/2. 

Next we consider the regime of strong non-Hermiticity. We choose Nf = 2 and detune the 
quark masses. The effect on the average sign is shown in Fig. 6. The sign problem is absent for 
mi = m2 and increases as |mi — mil increases. It again decreases with increasing v. 



5. Conclusions 

We have shown that a single random matrix theory describes two-color QCD at low density 
in the regime of weak-Hermiticity and at high density in the BCS superfluid phase, depending on 
the choice of the RMT parameter jx and on the rescaling factors in (2.2) and (2.5). These two 
regimes have very different symmetry breaking patterns. It would be interesting to investigate 
the applicability of random matrix theory in the region of intermediate densities, where intriguing 
phenomena such as a BEC-BCS crossover have been conjectured [16]. 

The analytical RMT results can be used to extract physical parameters such as A from lat- 
tice data. Two-color lattice simulations with adjoint staggered fermions, which are in the same 
symmetry class as continuum fundamental fermions [17], are currently underway to test our RMT 
predictions. 
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